'''
预制桩的下拉荷载(缺少)
'''

import math

d = float(input('桩直径(mm): '))
l = float(input('桩长(m): '))
l1 = float(input('桩进入第一层土厚度(m): '))
l2 = float(input('桩进入第二层土厚度(m): '))
l3 = float(input('桩进入第三层土厚度(m): '))
q_slk = float(input('请输入第一层土的极限侧阻力标准值(kPa): '))
lambda1 = float(input('请输入第一层土的抗拔系数: '))
q_s2k = float(input('请输入第二层土的极限侧阻力标准值(kPa): '))
lambda2 = float(input('请输入第二层土的抗拔系数: '))
q_s3k = float(input('请输入第三层土的极限侧阻力标准值(kPa): '))
lambda3 = float(input('请输入第三层土的抗拔系数: '))
qpk = float(input('请输入极限端阻力标准值(kPa): '))


# 单桩竖向极限抗压承载力标准值
if d<=0:
    print('请输入正确的桩直径')
elif 0<d<800:
    qsk = math.pi*d/1000*(q_slk*l1+q_s2k*l2+q_s3k*l3) 
    ap = math.pi*d**2/4/1000000
    q_uk = qsk + qpk*ap
    print(qsk,ap,q_uk)
else:
    phi_si1 = float(input('第一层土桩侧阻尺寸效应系数: '))
    phi_si2 = float(input('第二层土桩侧阻尺寸效应系数: '))
    phi_si3 = float(input('第三层土桩侧阻尺寸效应系数: '))
    phi_p = float(input('桩端阻尺寸效应系数: '))
    qsk = math.pi*d/1000*(phi_si1*q_slk*l1+phi_si2*q_s2k*l2+phi_si3*q_s3k*l3) 
    ap = math.pi*d**2/4/1000000
    q_uk = qsk + phi_p*qpk*ap

# 单桩抗拔极限承载力标准值
t_uk = math.pi*d/1000*(lambda1*q_slk*l1+lambda2*q_s2k*l2+lambda3*q_s3k*l3)

# 单桩竖向抗压承载力特征值
ra = q_uk/2

print('单桩竖向极限抗压承载力标准值为 %f kN' %(q_uk))
print('单桩抗拔极限承载力标准值为 %f kN' %(t_uk))
print('单桩竖向抗压承载力特征值为 %f kN' %(ra))